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ABSTRACT 

The conductance coefficients of disordered mesoscopic devices with n probes are in- 
vestigated within the noninteracting electron approximation at zero temperature. The 
probes are eliminated from the theoretical description at the expense of introducing non- 
local boundary conditions at the position of the contacts. Conductors with a large number 
of weak contacts are analyzed in detail. The ensemble-averaged conductance coefficients 
are in this case given by the so-called Hauser-Feshbach formula (gab) ~ PaPb where pa is 
the probability for emission of an electron into lead a. The conductance energy autocor- 
relation function is shown to deviate significantly from the conventional Lorentzian form 
and to have a negative tail. 
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1. Introduction 

At temperatures in the sub-Kelvin range, disordered metallic devices with a linear 
size L of about 1/xm or less display random but reproducible conductance fluctuations 
as a function of an externally applied magnetic field B [1]. It is now well understood 
[2] that these fluctuations, whose amplitude is universally of the order of e^//i, originate 
from the phase-coherent and diffusive quantum motion of electrons through the interior 
of the device, the "mesoscopic conductor". The sample-specific fluctuations with varying 
magnetic field are related, by an ergodicity argument [3], to statistical fluctuations across 
an ensemble of conductors characterized by the diffusion constant D. To estimate the 
correlation field one uses a scmiclassical path-sum argument and equates LP' /D with 
the time it takes in order for a typical diffusive electron path to acquire from the magnetic 
field Be an extra quantum mechanical phase equal to tt. 

The authors of Ref. [4] first pointed out that the value of Be is sensitive to the contact 
apertures joining the mesoscopic conductor to the current-voltage probes. They noticed 
that the mean time an electron dwells inside the conductor before escaping into the leads, 
is enhanced over Tdiff — L'^/D in the case of small apertures. Such an increase in dwelling 
time amounts to an increase in the area enclosed by a typical Feynman path, thereby 
leading to a reduction of Be- 

Ref. [4], while arriving at conclusions that are quite sound within the specific model 
considered, fell short of realizing the need for explicit introduction of a separate and inde- 
pendent time scale in addition to the diffusion time Tdiff. This step was made in Ref. [5], 
where the importance of the decay time Tdecay for the mesoscopic conductor problem was 
recognized. A lucid discussion of the physical meaning of Tdecay is given in Ref. [6]. The 
new time scale is set by the decay width F = ^/rdecay for electron emission into one of the 
leads. Ref. [6] quotes the formula 

F = A«iv, (1) 

where A is the level spacing of the mesoscopic conductor, is the number of open scat- 
tering channels at the Fermi energy Ep in the leads, and a is a dimensionless parameter 
measuring the quality of the contacts, (a ^ 1 for a good contact, and a ^ 1 for a weak 
contact.) When T ^ Ec = hD/L'^ or, equivalently, Tdecay <^ ''"diff, it is the Thouless 
energy E^ that defines the correlation energy of the problem and F becomes an irrelevant 
parameter. On the other hand, in the opposite limit F <S Ec, roles are interchanged and 
it is F, not Ec, that sets the relevant energy scale. Mesoscopic conductors satisfying the 
inequality T <^ E^ form the subject of the second half of this paper and will be called 
"decay- width dominated" for short. 

Formula (1) shows that there exist three independent ways of reaching the decay- width 
dominated regime: (i) reduce N by making the leads thinner; (ii) reduce a by degrading 
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the quality of the contacts; and (iii) if the conductor is three-dimensional (resp. quasi-one- 
dimensional), reduce A ~ L~'^ relative to ~ by increasing (resp. decreasing) L. 
The second option is pushed to the extreme by placing tunnelling barriers at the contacts, 
that is to say, potential barriers which electrons must tunnel through in order to reach the 
conductor from one of the leads or vice versa. 

The purpose of the present paper is twofold. First, wc wish to formulate the mesoscopic 
conductor problem in a way which, like the model of Ref. [5] , takes proper account of the 
existence of the decay width F, but differs from Ref. [5] in several respects: it does not 
exploit the connection with scattering theory but utilizes a more elementary expression for 
the conductance coefficients in terms of Green's functions; it uses the standard continuum 
Gaussian white noise model rather than Wegner's n-orbital model; and it eliminates the 
leads by a procedure known as R- matrix theory [7-10]. This procedure is particularly well 
suited for conductors with tunnelling barriers, where a treatment in terms of the scattering 
matrix seems awkward and may even be unfeasible. 

The second purpose is to communicate some analytical results for the conductance 
fluctuations of devices with tunnelling barriers and, more generally, of decay-width domi- 
nated conductors in the "locally- weak absorption limit" . By this term we mean the limit 
in which the electron emission rate per unit area of a contact becomes small and the total 
area of all contacts grows large, while their product is kept fixed. A special case of this 
limit are devices with a large number of weak contacts. To our knowledge, the limit of 
locally weak absorption has not been considered before in mesoscopic conductor physics 
(although its nuclear physics analog has received some attention [11]), perhaps because 
its experimental realizability is not clear. Nevertheless, we find it instructive to consider 
this limit since it (i) adds to our theoretical understanding of the phenomenon of univer- 
sal conductance fluctuations and (ii) stands out by some distinctive features concerning 
correlations. In fact, we will show that the functional form of the conductance energy 
autocorrelation function deviates markedly from the conventional Lorentzian one in this 
case. 

An outline of the contents of this paper is as follows. Sect. 2 defines the microscopic 
model and quotes the basic expression for the conductance coefficients gab {a,b = 1, n) 
in terms of advanced and retarded Green's functions at the Fermi energy. In Sect. 3, the 
problem of calculating gab is reformulated. The probes are eliminated and replaced by 
nonlocal boundary conditions on the surfaces separating the mesoscopic conductor from 
the leads. In Sect. 4, it is argued that for conductors with tunnelling barriers a local 
approximation for the integral kernel defining the boundary conditions may be used. This 
leads to the formulation of a phenomenological model, whose ensemble-averaged conduc- 
tance coefficients (gab) are calculated in Sect. 5 and shown to be given by the so-called 
Hauser-Feshbach formula known from the statistical theory of nuclear reactions. Sect. 6 
is concerned with the conductance covariance function {5gab{Ei)5gcd{E2)) ■ This function 
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turns negative as £'2! increases beyond F and satisfies the sum rule of having vanishing 
integral over E2 (or £^1). Our results are summarized in Sect. 7. 

I owe my understanding of many of the physical concepts and mathematical tools 
underlying this article to Hans A. WeidenmuUer. It is therefore appropriate that the 
article be dedicated to him on the occasion of his 60*^^ birthday and published in the 
present volume. 

2. The model and its conductance coefficients 

In this paper, three-dimensional disordered metallic devices of the generic type will 
be considered. Such devices consist of a mesoscopic conductor, roughly of linear size L, 
and of n leads joined to the conductor by the contacts. Both the conductor and its leads 
may have arbitrary geometrical shapes. 

Following common practice in mesoscopic physics, we assume that the statistical fea- 
tures of the conductance coefficients of such a device, at zero temperature, can be modelled 
by an ensemble of single-electron (or mean-field) Hamiltonians of the general form 

H=^{p-eAf + U + Um-(J+{p-eA)- (Uso x (j). (2) 

The functions U, and Vso are taken to be Gaussian white-noise potentials, with their 
strengths determined by the mean-free times for potential scattering, r, magnetic spin-flip 
scattering, r^, and spin-orbit scattering, T507 respectively [12]. The times r^, tso and 
the inverse cyclotron frequency (e/m x rotA)~^ are supposed to be long compared to r. To 
avoid complications that would otherwise appear in Sects. 5 and 6, we assume the contacts 
to be separated from each other by a distance of at least a few times the elastic mean-free 
path £ = vpT. 

The objects of the present investigation are the conductance coefficients Gab = e^/h x 
Qab (tt, b = 1, n), which determine the current response of the n-probe conductor to the 
applied electrostatic potentials. Our starting point is a formula for gab which originates 
from linear response theory and can be found in Ref. [13]; see also Ref. [14]. To write it 
down, the following definitions are needed. We take G'^ (resp. G~) to be the retarded 
(resp. advanced) one-particle Green's function of H at the Fermi energy Ep; formally: 

G^{x,y)= lim (EF±ie-Hy^{x,y). 

e— »0+ 

The limit e 0+ exists because the "openness" of the system makes the spectrum of 
H absolutely continuous. To keep the notation simple, we adopt the convention that the 
symbols x, y etc. comprise both position coordinates and spin projection. We introduce 
the conductivity tensor 

(^af3{x,y) = {vaG'^){x,y) {vpG~){y,x) + {G'^vp){x,y) {G~Va){y,x) 

+{vaG+v^){x,y) G~{y,x) + G+{x,y) {vpG'Vo^iy.x) (3a) 
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where Vq, is the a-component (a = 1, 2, 3) of the velocity operator 



dH 1 



(p - eA) + Uso X cr. 



(4) 



V = 



dp m 



For a G {1, n}, let Ca be any cross section of lead a. To account for spin, we introduce 
two identical copies of c^, one for each spin projection: c^d) and Ca(|), and we set Ca = 
Ca(T) U Ca(i). Eq. (75) of Ref. [13] can then be cast in the form (a ^ h) 



a,l3 

(Note that Ref. [13] ignores the spin degrees of freedom but the extension is straightfor- 
ward.) 

3. Elimination of the leads 

According to a standard assumption made in mesoscopic physics, the properties of 
the conductance coefficients gab are determined primarily by the phase-coherent diffusive 
motion of electrons inside the conductor and to a lesser extent by the details of what 
happens in the leads. In fact, one often takes the leads to be "ideal" or "clean" (although 
in reality the contacts and the leads are usually made from the very same dirty metal as is 
the mesoscopic conductor), i.e. one approximates the motion in the leads by free motion. 
Such an approximation is well justified if there is a clear geometrical division between the 
conductor and the leads, as effected for example by small contact apertures. In a situation 
where details of the motion in the leads do not matter, it makes sense to try and eliminate 
the leads from the theoretical formulation altogether. This can be done by a procedure 
known as R-matrix theory [9] in nuclear physics. (Note, however, that the roles played 
by the interior and the exterior of configuration space will be interchanged relative to the 
nuclear case.) It is applied in the present context as follows. 

Let a and b two fixed elements of the set {1, ...,n}. Evaluation of the conductance 
coefficient gat from Eq. (3) requires as input the Green's functions G'^{x,y) and G~{y,x), 
and their derivatives, for x E Ca and y E Cb- We will show how to construct G"*"; G~ can 
then be obtained from the relation G~{y,x) = G'^{x,y). 

We observe that the total configuration space, V, is partitioned by the surfaces Cc 
into the conductor space, Vq, and the lead spaces, Vc (c = l,...,n). Suppose now that 
we are to calculate G'^{x,y) for x E Ca and y E Vq. (We take y to Cb in the end.) To 




(3b) 
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begin, recall that G'^{-,y) = V'(-) is a solution of {Ep — H)'ip = on the region V — {y} 
with outgoing- wave boundary conditions. As as a first step, we take G~{z,x) to be the 
advanced Green's function of H (again evaluated at Ep) with support contained in Va 
and subject to the boundary condition G~{z, x) = for all z E Ca and x G Va- We write 
G~{-,x) = G~{-) for short. If V denotes the "covariant derivative" 

^ „ ie . im^^ 
V = V - -A + —Vso X cT, 
n n 

current conservation implies the validity of the relation 

div (g^ ■ V^j - VGx • V') = (5) 

on all of Va with the exception of the point x ^Va- The dots in Eq. (5) indicate summation 
over spin. Next, we integrate Eq. (5) over the region I4 — {x} and use Green's theorem to 
obtain ^ 

V'(^) = -|^/ {VaG^){z)'^{z)d^z, (6) 

where WDa/i — mVa with Va = Xla ^a"^" and is the vector normal to as before. 
Now we apply Va to both sides of (6) and take the point x to Ca- Using the relation 
Gx{z) = Ga{x, z) we then get 

{Vatp){x)= Ba{x,z)^{z)(fz {x E Ca) (7) 

where ^ 

Ba{x, ^) = y (^a^+Va) {x, z) {x, Z eCa). (8) 

Eq. (7) relates the normal component of the covariant derivative of i/j on Ca to the values 
and the derivatives of •0 on Ca- This relation can be regarded as a boundary condition 
satisfied by ijj- We have thus arrived at an exact reformulation of the original problem: 
instead of solving the Schrodinger equation for ip on the total space V — {y} (with outgoing- 
wave boundary conditions), we may solve {Ep — H)'ip = on the restricted space Vq — {v} 
supplemented with the boundary condition (7). The latter procedure produces all the 
information needed, viz. the values of = G'^{-,y) on Ca- 

Returning finally to complete notation and the specific task posed by Eqs. (3a-c), 
we observe that the Green's functions G^{p^q) for p^q & Vq are computed by solving 
{{^Ep — H)G'^){p, q) = 5{p — q) together with the boundary conditions 

{v^G+){u,q)= [ B^{u,z)G+{z,q)<fz = {G-Vc){q,u) (9a) 

{vcG-){u,q) = - [ B,{z,u)G-{z,q)£z = {G+Vc){q,u) (9b) 
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for all tt e Cc (c = 1, n). At the end of the computatfon we take p — > a; e and 
q ^ y & Cb- Let us now decompose the operators B,, as B^ = KeBc + iIm.Bc where both 
(ReBc){x,y) and {lmBc){x,y) satisfy the hermitecity relatfon 0{x,y) = 0{y,x). Eqs. (9) 
then permit us to recast (3b) in the form 

gah = n^ I (fx I d^y {{ReBa)G+){x,y){{ReBb)G-){y,x) (10a) 

JCa JCb 

where {c = a or c= b) 

{{ReBa)G^){x,y) = f {ReBc){x, z)G^{z,y)dh. (10b) 

JCc 

Of course, the usefulness of the above reformulation depends critically upon the pos- 
sibility of finding adequate and managable approximations to the nonlocal integral kernels 
Ba{x,y) (a = 1, ...,n). This possibility in turn is contingent upon our ability to controU 
the Green's functions see Eq. (8). 

In devising approximations to Ba, one should note that ImS^ corresponds, in a rough 
manner of speaking, to reflection at Ca- It is of minor importance and can often be 
neglected. On the other hand, ReBa causes "absorption" (taking the point of view of 
restricted configuration space Vq), that is to say, the loss of probability due to emission of 
electrons into lead a. This latter effect is crucial for determining the relevant time scales 
of the mesoscopic conductor, and it must therefore be modelled properly. 

By its definition in terms of through Eq. (8), Ba depends on the location of the 
surface Ca, and the difficulty of devising a good approximation varies with location. It 
is therefore very fortunate that current conservation leaves complete freedom in choosing 
the surfaces Ca, thereby permitting us to optimize the choice. How to make an optimal 
choice is rather obvious if there exists a clear geometrical division between the conductor 
and the leads, and if the leads can be taken to be clean: one will then place Ca right at 
the contact of lead a. In this case, an approximation which captures the essential features 
of Ba is 

Ba{x,y) = 7rJ2w^{x)Wc{^), (11) 

c 

where the sum runs over all open scattering channels at Ep in lead a, and the amplitudes 
Wa{x) are chosen phenomeno logically to fit the average rate of emission into the channels. 
This corresponds to the approximation used in Ref. [5] at the level of the S-matrix. Another 
case where Ba can be controlled with ease is analyzed in the sequel. 
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4. Reduced model for devices with tunnelling barriers 

From now on we will restrict ourselves to the case of n-probe conductors with weak 
contacts. In other words, emission into the leads will be supposed to be inhibited by ob- 
stacles of some kind. For definiteness, we imagine these obstacles in the form of tunnelling 
barriers and add a barrier potential Ub to the Hamiltonian H of Eq. (2). Let us denote 
by C'^ (resp. C") the surface separating the classically forbidden region under the barrier 
at contact a from the allowed region inside the conductor (resp. in lead a) . An optimal 
choice (in the sense of the previous section) is then to put Ca — C'^ (a = We 
adopt this choice. 

The simplification that occurs for devices of such kind is that one may use a local 
approximation to the kernels Ba- This is intuitively clear and can be quantified by the 
following semiclassical argument. 

Let a e {1, ...,n}. To construct Ba, we must first calculate G+(x,y) for two points x 
and y under the barrier and close to Ca, then apply the derivatives Va (see Eq. (8)) and 
finally take x and y to Ca- For simplicity, we will assume that the motion under the barrier 
is determined mostly by the interplay of Ub with the kinetic energy /2m and neglect 
the infiuence of all other terms in H. 

By applying the stationary-phase approximation [15] to Feynman's path integral for 
the propagator [16], we write G^{x,y) as a sum over classical paths with energy Ep 
connecting y with x. The contribution largest in magnitude comes from the path of shortest 
length, which, for Ub constant in the barrier region, would be a straight-line trajectory. 
This path, while generating the correct Green's function singularity as the distance |a; — y| 
goes to zero, is uninteresting here since it is an imaginary-time path without conjugate 
points and its contribution to G^ (resp. Ba) is purely real (resp. purely imaginary). 

An imaginary contribution to G^ arises from the path that takes off from y in the 
direction of lead a, bounces off the surface C" and then heads straight for a;. IHb denotes 
the thickness of the tunnelling barrier and Kp the average value of ^j2m(UB — Ef)/% 
along this path, the path's action is roughly given by 



There exist three more classical paths of a similar kind. These are the paths that are 
reflected at the surface Ca before heading for C'^ and/or before arriving at x. They carry 
opposite sign factors, and their combined contribution to the semiclassical approximation 
for G^{x, y) vanishes as x and/or y approaches Ca, as it should. 

In summary, by the above argument 



where da{-) is distance from Ca, and cq is a positive constant. Taking derivatives and 
sending x and y to Ca, we obtain 'RBBa{x,x) ~ exp(— 2^^^^), and we arrive at the 




l^G+{x,y) ~ -coda{x)da{y)e-^''^^'B+\^-y\'/^ 



8 



important conclusion that the range of ReBa{x, y) is of the order of ■\/TbT^- This range 
is to be compared to the elastic mean-free length, which is the smallest length scale 
of relevance for the diffusive motion inside the conductor. Inserting reasonable numbers 
{k~p^ > Ib ^ Inm and £ > lOnm), we see that i ^ ^/Ib/i^f- It is therefore justified, 
for the purpose of calculating the conductance coefficients, to suppress the length scale 
\/Ib/K'F and use a local approximation for Rei?a(a:, y). 

We are thus led to the following reduced model. We take reduced position space, Vb, 
to consist of the interior of the device only. OVq denotes the total boundary of Vq, and 

def 

C = U^^iCa (resp. dVo — C) its conducting (resp. insulating) part. We introduce a 
phenomenological function (3 : OVq — > R characterizing the penetrability of the barrier; 
P{x) ^ for X E C, and P{x) = for x e dVo — C. /? is assumed to be spin-independent 
(neglecting Zeeman splitting in the presence of a magnetic field, and magnetic spin-flip 
and spin-orbit scattering under the barrier). We impose the boundary condition * 

{vnip){x) = (3{x)ij{x) (aU X e dVo). (12) 

Here Vn is the component of v, Eq. (4), normal (in the spatial sense) to dVo- We take 
the Hamiltonian if to be a linear operator of the general form of Eq. (2), restricted to act 
on functions if) supported on Vq and subject to the boundary condition (12). We define 
G~^{x,y) to be the kernel of the operator {Ep — H)~^ with the boundary condition (12), 
and we set G~{x, y) = G~^{y, x). The conductance coeflicients are then calculated from 

9ab = h^ f d^x f d^y P{x)G+{x,y)P{y)G-{y,x). (13) 

JCa JCb 

This completes the definition of the reduced model. 

Of course, the Hamiltonian H, which was a hermitean operator in the original full 
space (including the leads, and with the usual square-integrability condition), becomes 
non-hermitean upon introduction of the boundary conditions (12), unless /? vanishes iden- 
tically. This non-hermitecity is an inevitable and, in fact, essential feature of any reduced 
description of the present type, and is particularly evident from the following alternative 
formulation of the reduced model. 

Let Sc be Dirac's ^-distribution with uniform support on C, i.e. 

/ Scfd''x= [ fd^x, 

JVo Jc 

* Although Dirichlet's boundary condition might seem more realistic for the insulating 
part of the surface, we feel that the internal consistency of the reduced model is enhanced 
by using boundary conditions of the form (12) everywhere. In any case, in the diffusive 
regime we intend to study, the conductance coefficients and their statistical properties are 
affected by the conductor-insulator boundary conditions only in a minor way. 
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and introduce an effective Hamiltonian 

H,f, = H--p6c. (14) 

Then , as defined above, is also the Green's function of the linear operator i/gflf acting 
on functions that satisfy the generalized Neumann boundary conditions 

To verify this statement, integrate {E — H^q)'^ = along an infinitesimal piece of curve 
intersecting C at a right angle, and use the boundary condition (15). From the continuity 
of ip across C, the boundary condition (12) is then recovered. The extra term —ihpSc/'^ 
in Hf.fl has an interpretation as an "absorptive contact potential", giving rise to the loss 
of probability which is caused by the escape of fiux through C. Note that the escape rate 
at a: e C, per unit area and spin projection, is given by I3{x)\'i(j{x)\'^ for a state with 
amplitude 'i/j{x). 

5. Average conductance coefficients of decay-width dominated devices in the 
locally-weak absorption limit 

Given an ensemble of one-particle Hamiltonians of the form specified in Sect. 2, we 
will now calculate for the reduced model of Sect. 4 the ensemble average (gab) (o 7^ b) for 
decay-width dominated conductors in the limit of locally weak absorption. 

Recall first the definition of decay-width dominated conductors by the inequality 
''diff <S '''decay where Tdiff = L"^ / D and, with vol = jy^ d^x, 

raLy = ^or' Ijd'x. 

Recall also the definition of the locally-weak absorption limit by /3 — > and d^x — > 00 
with j3d'^x kept fixed. A more quantitive condition for this limit to be attained is 

I3{x)f < / (5d^x for all xeC (16) 

JC 

where 1 = vpT is the elastic mean- free path for potential scattering. The inequality (16) is 
motivated by the observation that, since I is the lower cutoff length for difi^usive motion, 
d'^x should be compared to 

To compute {gah) from Eq. (13), we require the ensemble average 

P{x,y) = (|G+(x,y)|2) = (G+(x,y)G-(y,x)) 

for X E Ca and y E C},. We will now calculate this quantity in the prescribed limit. To 
begin, let x and y be any two points in Vq. The inequality Tdifj <^ ^decay implies that an 
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electron moving diffusively inside the mesoscopic conductor, traverses the distance L many- 
times before being absorbed at one of the contacts (that is to say, before being emitted 
into one of the leads). The density \G~^{x, y)\'^ of a stationary state sustained by a source 
at y therefore becomes independent of x upon ensemble averaging, with two exceptions. 
The first occurs whenever x is within a distance of £ or less from one of the contacts, where 
wave functions, and consequently P{x, y) as well, are depressed by absorption in general. 
However, the condition (16) for absorption to be locally weak precisely means that this 
depression can be neglected under the present circumstances. The second exception occurs 
when X lies within a distance of I or less from the source point y. In this case, the value 
of P{x, y) is modified by coherent backs cattering, as is well known from weak localization 
physics [17]. Summarizing this paragraph, we write 

Pix,y) = Pd + Pcix,y) (17) 

where P^ is the constant part of P and Pc is the short-ranged backscattering correction. 
The sign of Pc{x,x) depends on the relative strengths of potential scattering, magnetic 
spin-flip scattering, spin-orbit scattering and the magnetic field [12], but its magnitude 
is always comparable to Pd- If the arguments presented in support of Eq. (17) seem too 
heuristic, a technical derivation is given in Appendix A. 

To proceed, recall that we assume the contacts a and b to be separated from each 
other by a distance of at least a few times i. By Eq. (13) we then need to know P^ but 
not Pc{x,y). Multiplying both sides of Eq. (17) with P{x), integrating over J^d^x, and 
dividing by (3d'^x, we have 

Pd - / d^x Pix){Pix,y)-Pc{x,y))/ [ 0x. (18) 

JC JC 

The contribution to the right-hand side of (18) from Pc is negligible under the condition 
(16), since 

/ d\ Pix)Pcix, y) ~ fmPciv, y) « Pd I ^d\. 
Jc Jc 

The remaining term can be calculated by using the identity 

in f d''yG+{x,y)/3{y)G-{y,x)^{G--G+){x,x), (19) 
Jc 

which is a consequence of the boundary condition (12) and of current conservation. Taking 
the ensemble average on both sides and using 

{{G- -G+){x,x)) =27rw, 

where v (independent of x) is the local density of states, we obtain 
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We insert this relation into the ensemble-averaged version of Eq. (13). Upon introduction 
of the quantities 



n 



a= — 

vol 



A-i = zyxvoi, r = ^r, 



a = 1, ...,n), 



0=1 



where A -'^ is the total density of states (counting spin), T the total decay width, and Ta 
the partial decay width for emission into lead a, the expression for (gab) takes the form 



This is the analog of what is called the "Hauser-Feshbach formula" in the statistical theory 
of nuclear reactions [18]. Note the interpretation of the ratio Ta/T = (3cPx/ J^f3(Px 
as the probability for emission of an electron from the mesoscopic conductor into lead a. 
Note also that in comparison with the Thouless formula [19], the Thouless energy has been 
replaced by 27rTaTi,/T. 

Let us mention in passing that the validity of the Hauser-Feshbach formula is not 
confined to the locally-weak absorption limit. As a matter of fact, all that is needed in 
order for (gab) to have the factorized form of Eq. (20), is the long dwelling time of a 
decay-width dominated conductor, causing the processes of entry from lead a (or b) and 
emission into lead b (or a) to be uncorrelated. Thus, formula (20) remains valid for strong 
absorption if the partial decay widths F^ are replaced by more complicated, nonlinear 
expressions in We will not elaborate upon this point here. 

6. Conductance covariance function of decay-width dominated devices in the 
locally-weak absorption limit 

Considering devices of the same special kind as before, we will now calculate the 
correlation function 



for a 7^ 6 and c ^ d. This will be done for the reduced model of Sect. 4, taking /3 to 
be energy-independent. * Eventually, we will specialize to the three universality classes 
[20,17,12] that are known to exist for the mesoscopic conductor problem. Recall [21] that 
these are denoted by I (orthogonal class: potential scattering), Ila (unitary class: magnetic 
field), lib (unitary class: magnetic spin- flip scattering) and III (symplectic class: spin-orbit 
scattering). 

* The justification for neglecting variations of gab due to changes in the penetrability 
coefficient (3 - which are rather strong in a microscopic model with tunnelling barriers - is 
that such variations do not concern us here. What is measured in experiments with variable 
magnetic fields are the disorder- induced statistical vairiaiions of gab, and it is therefore the 
latter that we wish to calculate. 




(20) 



{gab{EF + dE)gUEF)) 
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With all the preparations made in Sect. 5, we can now be rather brief. Proceeding as 
before, we are led to consider the ensemble average 

{G^{xa, Xb; Ei)G~{xb, Xa, Ei)G+{xc, xa] E2)G~{xd, Xc, E2)) (21) 

for a^o e Ca, Xd e Crf. Here it is necessary to distinguish cases. Let us first assume that 
the indices a, 6, c and d are all mutually different and, as before, contacts are separated 
from each other by a distance of at least a few times the elastic mean-free path I. Then the 
argument of Sect. 5 goes through without change and expression (21) can be approximated 
with negligible error by a constant, independent of Xa-, Xd- Introducing the conductance 
covariance function 



Cab,cd{Ei,E2) = {gab{Ei)g,d{E2)) - {gab{Ei)) {g,d{E2)) 
and using again Eq. (19), we find 

Cab,cd{Ei, E2) = {gab{Ei)gcd{E2))R{Eu ^2), (22a) 

where 



R(F {l^G+{E^)xlmG+{E2)) 

^ {ImG+{Ei)){lmG+{E2)) ' 



(22b) 



and G'^{E) = G (E) is the trace of the Green's function of the Hamiltonian H with the 
boundary condition (12). 

If, on the other hand, at least two indices in (21) coincide, or if contacts are sepa- 
rated by a distance less than £, a complication arises. Consider for definiteness the case 
a = d. There now arise contributions to the double integral d^Xa d'^Xd from the 
region \xa — Xd\ < In this region, the value of the otherwise constant expression (21) 
is modified by cross contractions between G^{xa-, Xb] Ei) and G~{xd, Xc] E2). In diagram- 
matic language, these give rise to impurity ladders of the same "cooperon" type that cause 
the backscattering correction discussed earlier. (For a — c diffuson ladders appear, too.) 
Because of the condition \xa — Xd\ < the additional terms lead to contributions carrying 
the factor /3'^d'^x, which is to be compared to the regular contribution, Eq. (22a), 

carrying the factor (J^ fid'^x)'^. Thus the terms missing in Eq. (22a) are small if 

/ r \ 2 

32-2, 



/3^d^x<^[ I3d^x] (a= 1, ...,ri). (23) 

Note that this condition is much more restrictive than (16). In summary, in order for 
Eq. (22) to be correct, we must either suppose that the indices a, b, c and d all differ from 
one another (and \xi — Xj\ > £ for all Xi e Cj, xj e Cj (i ^ j)) or else impose the stronger 
condition (23). 
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Eq. (22) reduces the problem of calculating Cab,cd{Ei, E2) to a problem in "level 
statistics" for an open system. Note that in contrast to Ref. [19], where a somewhat 
similar relation appears, Cab,cd is here expressed directly in terms of R{Ei, E2) and not 
as a double integral over energy. The possibility to do so is offered by the limit under 
consideration. 

The calculation of (22b) is still not easy in general but in the locally-weak absorption 
limit and for each of the universality classes, it reduces to a problem whose solution is 
known. For the following it will be convenient to represent G^{E) = G^{E) alternatively 
by G'^{E) = iT{E—Hf.Q)~^ with i^efr = H—ihpSc/'^ and the Neumann boundary condition 



Consider the matrix element M^i of h/SSc between two eigenstates ipk and ipi of H: 



Jc 

The phases of t/jk and t/ji are randomized over a correlation length of the order of the elastic 
mean-free path £. We therefore expect 



These estimates mean that hl35c acts effectively as a multiple of the unit operator under the 
condition (16). In other words, for the purpose of calculating R{Ei, E2) we are permitted to 
approximate G'^{E) in Eq. (22b) by ii{E-\-iT /2 — H)~^. To avoid possible misconceptions, 
we stress that the approximation %f35c ~ F x 1 must not be over interpreted. For the case 
F ^ A, which is included in the limit where our considerations apply, the operator fi(35c 
docs cause appreciable mixing between states close in energy. However, such mixing does 
not modify the correlation properties of the spectrum. In particular, it does not reduce the 
spectral rigidity. The approximation is therefore justified when used for the purpose of 
calculating correlation functions such as R{Ei,E2). 

For the special case F = 0, the two-level correlation function resulting from the sub- 
stitution G~^{E) — > tr{E + ir/2 — H)~^ in Eq. (22b) has been calculated exactly for each 
universality class by Efetov [22] in his work on the level statistics of small metallic particles. 
It is straightforward to extend Efetov's calculation to F 7^ 0. The result for R{Ei,E2), 
obtained by performing this analytic continuation, is quoted in Appendix B. 

Here we refrain from considering the general case but specialize to F A, the case of 
an "open" system, which is also accessible (without Efetov's results) via standard diagram- 
matic perturbation theory. Using either the exact result of Appendix B, or more simply 



(15). 





and, by the law of large numbers. 
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the perturbative cooperon-difFuson expansion - see Eq. (33) of Ref. [19] - one obtains 



E2) = Re (El -E2 + iV)-^ , (24) 

where the distinction between universahty classes is contained in the coefficients ci = 16, 
ciia = 8, Cm = 4 and cub = 2. Insertion of (24) into (22a) and use of Eq. (20) gives 

1 — (E — E \^ /r^ 

Cab,cd{El, E2) = CkPaPbPcPd (^l^(^Ei-E2)^/T^)^ ^^^^ 

where pa = Ta/T, pt = Tb/F etc. 

The universal numbers Ck [k — I,..., Ill) reflect the influence of symmetries on the 
conductance fluctuations. Breaking of time-reversal invariance (I Ila or III lib) 
causes a reduction by a factor of 2, and breaking of spin-rotation invariance (I 111 or 
Ila — > lib) causes a reduction by a factor of 4. These reduction factors are not new but 
were already observed in [19]. 

What is more striking about Eq. (25) is the non-Lorentzian dependence of Cab,cd on 
(£"1 — E2)/r and the validity of the sum rule 

/ CabME,E')dE' ^0, (26) 

which is seen most easily from Eq. (24) and Cauchy's theorem. We emphasize that the 
result (26) does not hold in general but is a special feature of the locally-weak absorption 
limit. (It is true, however, for any value of the decay width T of decay- width dominated 
conductors; see Appendix B.) To understand this feature, we recall that the rigidity of the 
energy spectrum of an isolated system is preserved for an open (but decay-width domi- 
nated) system satisfying the condition (16): energy levels are, roughly speaking, shifted by 
a common amount E ^ E — iV /2. The negative tail, or "correlation hole", of R{Ei, E2) 
for \Ei — E2\ > r is a direct consequence of this rigidity. 

The conductance fluctuations result from Eq. (25) by setting Ei = E2: 

(Sgab^gcd) = CkPaPbPcPd- 

These fluctuations are universal in the sense that they are sensitive - leaving aside the 
dependence on symmetry - only to the geometry of the n-probe conductor as expressed by 
the coefficients (^ = 1, ?i). However, this is not a universality in the strictest sense of 
the word since Pa — Ta/ ^^=1 T/c does become smaller as the total number n of contacts 
increases. Strict universality is recovered by dividing the contacts into two groups, say 
"left" (L) and "right" (R), and summing over a, c e L (a 7^ c) and b,dER{b^d). 
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7. Summary 

This paper was in two parts. In the first part, we addressed the practical question of 
how to calculate the conductance coefficients gab (a, b = 1, ...,n) from the linear-response 
formula (3a-c). Supposing that the statistical properties of gab are determined mostly by 
the diffusive motion of electrons inside the mesoscopic conductor, it is natural to try and 
eliminate the leads from the theoretical description in favour of suitable boundary condi- 
tions imposed on the contact surfaces Ca {a — 1, ...,n). This idea is usually implemented in 
a very schematic fashion. Refs. [2-4] invoke the condition Pd = (resp. V^-Pd = 0) on the 
conducting (resp. insulating ) part of the surface, at the level of contructing the diflFuson 
propagator Pd. Such a condition, while reasonable for thick leads and good contacts, is 
too crude for thin leads and/or weak contacts. (The authors of Ref. [4] are of course aware 
of this limitation. This is why for the thin-lead case they adopt the trick of making a 
somewhat arbitrary division of the leads into clean and dirty regions.) 

Inspection of Eqs. (3a-c) shows that what we need to do is to calculate the covariant 
derivatives {vaG^){x,y) for x e (a = l,...,n). In Sect. 3 we restated the problem of 
making this calculation on the total space (including the leads) as an exactly equivalent 
but reduced problem formulated on the bounded domain of the mesoscopic conductor. 
The reduced formulation involves a set of integral kernels Ba, see Eq. (8), which define 
boundary conditions expressing the covariant derivative of normal to Ca in terms of the 
values of on Ca (a = 1, ...,?!). These nonlocal boundary conditions cause the reduced 
(or effective) Hamiltonian to be non-hermitean and energy- dependent. Ba is determined 
exclusively by the properties of electronic motion in the leads and is therefore easy to 
controU when the leads are clean; see Eq. (11). (We did note translate the boundary 
conditions (9) with Ba given by (11) into a corresponding boundary condition for the 
diffuson propagator Pd- How this is done is described at length for a different microscopic 
model in Ref. [5]; see also Sect. 7 of Ref. [23] for some useful mathematical details.) 

Another example where Ba can be controlled are conductors with tunnelling barriers. 
The semiclassical argument given in Sect. 4 shows that a local approximation to Ba may 
be used in this case. This argument led to the formulation of a phenomenological model, 
whose diffuson propagator Pd satisfies 



with the boundary condition —DVnPd — PPd on the surface C of the conductor. 

In the second part of the paper, the phenomenological model of Sect. 4 was analyzed 
for decay-width dominated conductors in the locally-weak absorption limit. The technical 
simplification occurring for decay-width dominated conductors is that the spectral expan- 
sion of {—?iDV'^)~^ is dominated by the smallest eigenvalue, Aq, of —TiDV^ (with the 
specified boundary condition) in this case. Aq and the next eigenvalue, Ai, are given by 



no {V^Pd){x,y) = 27riy 5{x - y) 




V + const X Ec {Ec = hD/L^). 
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Note that for a three-dimensional mesoscopic conductor the hmit Aq <^ Ai is approached 
by keeping the conducting part of the surface constant in size and making L larger. 

The focally-weak absorption hmit is defined technicaUy by the condition (16). It 
guarantees the correlation properties of the eigenvalues of the effective Hamiltonian H^.^ = 
H — ihpSc/'i to be the same as those of H — iT/2. In other words, the eigenvalues of the 
Hamiltonian H of the isolated system are simply shifted, roughly speaking, by a constant 
amount —ir/2. 

By expressing the conductance fluctuations as a density-of-states correlation function 
in Eq. (22) - which is rigorously justified in the prescribed limit and for a, b, c, d all 
mutually different - we obtained 



(l + (Ei-E2)Vr2)2- 

The interesting feature of this expression is that it turns negative at \Ei — > T. The 
negative tail can be understood as being a consequence of the spectral rigidity of H. 

Finally, let us translate the correlation energy F into the corresponding magnetic-field 
correlation scale B^. measurable in experiments on a single sample. To do that, we will 
follow a procedure described in Ref . [4] . A difference in magnetic field B — > B + couples 
into the equation for the diffuson propagator by the substitution V — > V — ie5A/h where 
5B = TotSA. We need to calculate the correction to Aq in lowest order of 5B. Perturbation 
theory gives Aq — ^ Aq + const x Ec{6(p/(po)'^ where d(p/(pQ is the change in magnetic flux 
through the conductor, measured in units of the fiiix quantum 0o — h/e. By equating the 
correction to F, we obtain 

5eL2~(/>ov^f7^. 

Thus, in comparison with the case of thick-lead conductors with good contacts. Be is 
reduced by a factor of ^JV jE^. 
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Appendix A 

To establish Eq. (17) on a rigorous level, we find it convenient to make use of the 
elegant and powerful technique of Efetov [21], who has shown that P{x,y) = {\G'^{x,y)\'^) 
can - in the diffusive regime - be regarded as the propagator of certain Q-matrix superfields 
with "free energy" functional It is not difficult to adapt Efetov's treatment to the 

reduced model with boundary conditions (12) or the alternative formulation in terms of 
the effective Hamiltonian (14). For simplicity, we assume that electrons are spinless and 
subject to potential scattering only. In this case, one finds that Eq. (3.53) of Ref. [21] is 
replaced by 

^[Q] = stvQ'{x)d'x+^ {sir\n{EF _ |!. + l^Q + |/?A5c)) {x,x)d^x, 

and P{x, y) is a sum of diffuson [d) and cooperon (c) contributions, P = + Pc-, where 

Pd{x,y) = Jt^Q Gri^^^-^Q)Gd^iy^y^Q)^~^^'^\ (A-la) 
Pc{x,y) - Ji^Q gt~i^^y,Q)Gc^iy,^-,Q)e-^^'^\ (A.lb) 



and 



g{x, y;Q)^(EF-^+'^Q+ '-^(3A6c) {x, y). 



The notations used are either self-explanatory or those of Ref. [21]. In the next step, one 
makes the usual saddle-point approximation followed by a gradient expansion [21], valid 
in the metallic regime kpi ^ 1. Retaining in J-' only terms up to linear order in /3, one 
obtains 

J^[Q] c^^(^D J str{VQfd^x + 2 J PstrAQd^x^ (A.2) 
where Q satisfies the nonlinear constraint = 1- Eq. (A.2) omits such terms as 



const X {hnuff / sir {AQfd'^x. 
'c 



These are negligible under the locally- weak absorption condition (16). With the same 
degree of accuracy, one may use the approximation ^(x, x] Q) ~ —mvQ^x)^ which leads to 

Pa{x,y) = -tt'i^^J VQ Q+- {x)Q^+ {y)e-^^Ql 

To analyze this expression further, one observes that the condition Tdifj <^ Tdecay for a 
conductor to be decay-width dominated can be written 

hD h f 
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A glance at Eq. (A. 2) then shows that the dominant contribution to Pd{x,y) comes from 
spatially constant Q-fields, spatially varying Q-fields being separated from the constant 
ones by a large gap in free energy. This proves the independence of P(i{x,y) of x and y. 

The other properties used in Sect. 5 are Pc{x, x) ~ Pd and the short range of Pc{x, y). 
These follow from Eq. (A. lb) and 

g{x, y; Q)^(^Ef-^ + (x, y) f{x - y)Q{{x + y)/2) 

where /(O) = —iniy, and / decays over a length scale of the order of the elastic mean- free 
path. 

Finally, we observe that the various perturbations in H, Eq. (2), which break time- 
reversal and spin-rotation symmetry, do not affect the diffuson degrees of freedom, whose 
free energy is always given by a functional of the form (A. 2). Therefore, P^ is constant in 
general. 

Appendix B 

With T = h J^/3d'^x/vo\ the decay width, let pr{E), 

p^{E) = --Im tT{E + iT/2 -H)-^, 

TT 

be the total density of states (counting spin) of the open system. We will write down exact 
expressions for the two-level correlation function 

for all universality classes, obtained by analytic continuation of Efetov's results [22]. 
We introduce the functions (C — > C) 

= (B.la) 

\7r OZ ITTZ J Jq t 

and, with = {py{Ef))-, we set 

Fi{z^) = f{z/2) + g{z/2), 

Fu,izA) = f{z/2), FuhizA) = f{z), (B.2) 
Fui{zA) = f{z) + h{z). 
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R is then given by 

Rk{E^,E2; r) = ^{Fk{E^ -E2 + iT) + F^iE^ -E^ + iV)) (B.3) 

in all cases {k = I, Ila, lib, III). 

Note that the function /, Eq. (B.la), is holomorphic for ^ e C — {0} and vanishes as 
z~'^ at infinity in the upper half of the complex plane. Hence, by Cauchy's theorem, 

f{x + i^)dx = for 7 > 0. 

The same statements apply to the functions g and /i, Eqs. (B.2b) and (B.2c), which leads 
to the sum rule 

/ R{E,E';r)dE' = 0. 
Jn 

The result (24) for £^2; T) quoted in Sect. 6 follows from the asymptotic expansions 

*'(n) = ^(l + ;^ + 0(7-)), 

M.7) = ^(l-|^+0(7-)), 

for 7^1. Note that Eq. (B.4a) shows the corrections to Rk, Eq. (24), for k = Ila 
and lib to be exponentially small and therefore not calculable by perturbation expansion 
in [El — E2 + ir)~^. Furthermore, from the exact expressions (B.1-3) one can obtain 
the conductance fluctuations {SgabSgcd) for small values of F, which are not accessible 
by perturbation theory either. Both f{i^) and g{ij) behave for 7 — > as {7r^)~^, while 
h{i^) ~ |7r7ln7r7. Insertion of these limiting forms into Eqs. (B.2), (B.3) and (22a) gives 

i^gab^gcd) ^ const X PaPbPcPd X F/A. 

This shows that the conductance fluctuations go to zero for F — > 0, which is the limit of an 
isolated system. In other words, the behavior of (SgabSgcd) is dominated by the decrease 
of {gab){gcd) ~ F^/A^ for F — > 0, overpowering the increase in relative fluctuations due to 
the appearance of isolated resonances. This conclusion is contrary to the one reached in 
Ref. [4] . We attribute the discrepancy to overextension, in Ref. [4] , of perturbation theory 
to a regime where it does not apply. 
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(B.4a) 
(B.4b) 
(B.4c) 
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